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Abstract 



^ ■ Fractional dynamics of relativistic particle is discussed. Derivatives of fractional orders 

Q^l with respect to proper time describe long-term memory effects that correspond to intrinsic 

Q . dissipative processes. Relativistic particle subjected to a non-potential four-force is consid- 

^ ' ered as a nonholonomic system. The nonholonomic constraint in four- dimensional space- 

>-C , time represents the relativistic invariance by the equation for four-velocity u^j,u^ -|- c^ = 0, 

where c is a speed of light in vacuum. In the general case, the fractional dynamics of rela- 



tivistic particle is described as non-Hamiltonian and dissipative. Conditions for fractional 
relativistic particle to be a Hamiltonian system are considered. 
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1 Introduction 

The derivatives of non-integer orders are a natural generalization of the ordinary differentiation 
of integer order. Fractional differentiation with respect to time is characterized by long-term 
memory effects. The theory of derivatives and integrals of non-integer order goes back to 
Leibniz, Liouville, Riemann, Grunwald, and Letnikov PQ|2]. The interest in fractional equations 
[21 [3] has been growing continually during the last few years because of numerous apphcations 
in recent studies in mechanics and physics (for example, see books [H [5l [6] and references 
therein) . 

In this paper, we discuss fractional dynamics of relativistic particles that are described as 
nonholonomic systems in four-dimensional space-time. It is well known that components of 
the four- velocity u^ = dx^/dr (yU = 1, 2, 3, 4 and r is a proper time) are not independent. The 
components of the four- velocity are connected by the equation u^u^ + c^ = 0, where c is a speed 
of light in vacuum. This equation allows us to consider the relativistic particle as a system with 
constraint in four-dimensional space-time. This constraint is nonlinear nonholonomic (nonin- 
tegrable) constraint. As a result, the relativistic invariance for point particles is represented by 
a nonholonomic constraint [3 El [9] . 

Note that only mechanics of relativistic particles can be considered as a mechanics with 
nonholonomic constraint. The relativistic invariance in the field theory cannot be represented as 
a nonholonomic constraint. At the same time, nonholonomic constraints can be used in the field 
theory. For example, the higher spin fields are connected with nonholonomic constraints [lO] and 
the gauge fixing conditions for of non-abelian gauge fields can be described as nonholonomic 
constraints [11]. The Euler-Lagrange and Hamilton equations for nonholonomic systems in 
classical field theory are suggested in [T2] . 

In the framework of the fractional dynamics, we consider a relativistic particle subjected 
to a general four-force. In the general case, the four-force is non-potential, and the relativistic 
particle is a non-Hamiltonian system in four-dimensional pseudo-Euclidean space-time. We 
consider fractional dynamics of non-Hamiltonian and dissipative systems in relativistic the- 
ory. The fractional equations of motion describes power-law memory effects that correspond 



to intrinsic dissipative processes in the relativistic systems. Note that relativistic particle with 
dissipation is discussed in fi3\ HI]. In Refs. [131 E]) the Lagrangian and Hamiltonian func- 
tions for one- dimensional relativistic particles with linear dissipation are suggested. In general, 
non-Hamiltonian and dissipative n-dimensional systems with n > 1 cannot be described by 
Hamiltonian or Lagrangian since the Helmholtz's conditions for these systems are not satis- 
fied [15]. In this paper, we consider fractional dynamics of relativistic particles as motions of 
four- dimensional non-Hamiltonian and dissipative systems. 

In Sect. 2, the nonholonomic constraint in four- dimensional space-time for relativistic par- 
ticle and some notations are considered. In Sect. 3, we discuss the fractional equations of 
motion for relativistic particle that is considered as a nonholonomic system. In Sect. 4, we 
discuss the d'Alembert-Lagrange principle for fractional equations of relativistic particle that 
is considered as a nonholonomic system. We prove that fractional equations for relativistic 
systems with nonholonomic constraint are represented as fractional equations for holonomic 
systems. In Sect. 5, the conditions for fractional relativistic particle to be a Hamiltonian or 
non-dissipative system are considered. Finally, a short conclusion is given in Sect. 6. 

2 Nonholonomic Constraint 

We consider a four-dimensional pseudo-Euclidean space-time of points with coordinates x^: 
x^ = X, x'^ = y, x^ = z, x'^ = ct. The point coordinates in the four- dimensional space-time 
can be considered as components radius four-vector of the point particle, R = (x^, x^, x^, x'^) = 
(x, y, z, ct). The square of the elementary radius four- vector in the four-dimensional space-time 
is defined by {dRY = rj^j^^dx^dx'^ . Here and later we mean the sum on the repeated indices /x 
and V from 1 to 4. The coefficients rj^y define a metric of pseudo-Euclidean space-time. This 
metric is a diagonal tensor such that r/n = 7722 = "/^ss = 1 and r/44 = — 1. Note that x^ is not 
equal to x'^, since x^ = tj^^x" and xi = X"*^, X2 = x^, X3 = x^, and X4 = — x'^. 

Assume that we have two radius four- vectors R and R' with coordinates x^ and x'^ of two 
reference frames to describe a relativistic particle. If the coordinate transformation x'^ = a^^x", 



where a^ are constant values, satisfies the invariant condition: 

{dE!f = {dRf ■ T^^M^dx"" = 7]^^dx''dx^, (1) 

then this transformation is a Lorenz transformation. The coordinates of the radius four-vector 
in the proper reference frame are Rq = (0,0,0,cr), where r is a proper time. Condition (JT]) 
leads us to the relation 

(dR)'^ = {dRof : r]^^dx''dx'' = -c^dr'^. (2) 

Using the definition of three-velocity v'^ = dx'^/dt, k = 1,2, 3, we get 

dt = ^dT, ^= (l-vycY'^\ (3) 

Four-velocity of the point particle is defined as a derivative of the radius four-vector with 

respect to proper time: 

- dR „ dxt" 
V = ^ : u^ 



dr dr 

dr 1^ ■: ■" ^: ^: '^y ""^ '^ dr 



The components of the four-velocity V are u'^ = ^ = •yv'^, k = 1,2,3, and u^ = ^ = cy. 



Note that rest particles {v = 0) have u^ = c. 
Equation ([2]) leads to the relation 

'dR\ _ fdRo\ dxi'dx" _ ^ 

dr j \ dr j ^^ dr dr 

which means that square of the four- velocity is a constant value: V"^ = —(? . Therefore we have 
the constraint equation 

^^.u^v" + c' = 0. (4) 

As a result, a relativistic particle in the covariant formulation of relativistic mechanics is a 
system with the nonholonomic constraint. The constraint (jlj) is nonholonomic since it depends 
of velocity. Relativistic mechanics can be considered as nonholonomic mechanics in the four- 
dimensional space. 

It is known that constraints in mechanics are some simplifications of real particle interac- 
tions. (Note that this statement is not correct in the field theory.) For example, if we consider 



the pendulum then we usually neglect of the forces of thread deformation. We also neglect of 
an interaction for constraint (|4]), which defines the relativistic invariance. If we use the non- 
holonomic constraint (j4]), then we neglect of a gravity interaction between particles. Let us 
consider the deformation of equations ([2]), ([3]) and (jlj) in general theory of relativity [16]. In 
the approximation of weak gravity fields, we have 

{dRf = ri^^dx'^dx'' - 2ipdt^ = -c^dr^, 

where 

dt = j'dT, 7' = (^1 + -^ - _ j , 

and ip is a. classical (Newtonian) gravity potential. As a result, we have 

Therefore nonholonomic constraint 1^, which defines the relativistic invariance, is connected 
with the neglect of the gravity interaction, ip = 0, (in general theory of relativity). 

3 Fractional equations of motion of relativistic particle 

Let rriQ be a rest mass of a point relativistic particle. The four-momentum of the particle is 
defined by P = moV . The components of the four-momentum are p^ = ttiqU^. Equation (jlj) 
gives 

Vi.^p'^p" + mlc^ = 0, (5) 

In relativistic mechanics the Newtonian equations are replaced by some generalization, which 
is invariant under the Lorenz transformations [171 [16] . The Newtonian equations are satisfied 
in the proper reference frame. The four-vector analog of the Newtonian equations is 

dP ^ ^ ^ 

— = J^(r,R,P). (6) 

This equation is postulated as a main equation of relativistic dynamics. Equation ([6]) describes 
a relativistic particle subjected to a four-force ^ = F{t, R, P). Equation (Q must be considered 



with condition (E]). As a result, we have the equations 

-— = — p^ -^ = J-^(r, x,p), r]^,p^p'' + mlc^ = 0. (7) 

dT mo dr 

If dmo/dr = 0, then (JTj) give 

moD^x'' = J^''(t,x,p). (8) 

where ri^yD].x^D].x^ = —(? . These equations of motion can be generahzed for fractional dy- 
namics to take into account a power-law memory. We consider a generalization of ([H]) in the 
form of the fractional differential equations 



mo 



^D^x^ = J-^(r,x,p), (1< a < 2) (9) 



involving the Caputo fractional derivative f -D". The left-sided Caputo fractional derivative [2] 
of order a > is defined by 

° ^-"^ " V{n - a) io (r - rO'^'^+i " °^ ^ ' ^ ^ 

where ra — 1 < a < n, Z^" = d^/dr^, and o-^r is the left-sided Riemann-Liouville fractional 
integral 

oOM = ^r/^^^, (r>0). (11) 

r(«) io (^ - ^ ) 

Fractional derivative with respect to proper time describes a power-law memory effects that 
correspond to intrinsic dissipative processes. 

Using p^ = rrioDlx'^, equations ([9]) can be rewritten in the form 

Dlx^" = — p^ (12) 

mo 

^DrV = ^^{r,x,p), (l<a<2), (13) 

Fractional integration of flT3|) of order a — 1 gives 

o/r^^z^rv = oir'j"'ir,x,p). (m) 

Using the fundamental theorem of fractional calculus [18] 

oir' oDrV = P^{r) - p(0), (0 < 1 - a < 1), 



we obtain 

p'^{T)=p'^iO)+oir'J"'ir,x,p) (15) 

Differentiation of (fT5|) gives 

Dy = o/^^"-F'^(r,x,p), (0 < 2 - a < 1), (16) 

where o-D^~" is the left-sided Riemann-Liouville fractional derivative defined by 

oD^x^ = D': o/r°a;^ = ft V"^ / 7 ^^ ZT> (^ - K a < n). (17) 

As a result, equation ([9]) is equivalent to the fractional equations 

mo 
Z^y = o/^'""^'^(r,x,p), (1< a < 2). (19) 

These equations describe fractional dynamics of relativistic particle. Fractional differentiation 
with respect to proper time is characterized by long-term memory effects that correspond to 
intrinsic dissipative processes in the relativistic systems. 

4 d'Alembert-Lagrange principle for fractional relativis- 
tic dynamics 

It is known that the general principle, which allows us to derive equations of motion with 
holonomic and nonholonomic constraints, is the d'Alembert-Lagrange principle. For equations 
(IT8|) and (IT9|) this principle leads to the variation equation 

(^ - ^Dl-^F^{T,x,p))r^,M'' = 0. (20) 

Multiplying (iT9l) on the variation fe^ = rj^^dx'^ and summing over yU we obtain this variational 
equation. 

The variations of coordinates 5x^, /i = 1,...,4 are defined by the relation of the ideal 
constraint 

7^^5x^ = 0, (21) 



where TZ^ are components of the constraint force vector. The four- vector 7?.^ can be considered 
as a contribution of the reaction associated with the constraint to the four- force o-D^^"J-'^(r, x,p). 
Because a reaction force does no work in a virtual movement that is consistent with the corre- 
sponding kinematical restriction, we conclude that TZ^ must be perpendicular to any dx'^ that 
satisfies the constraint equation. Thus, if 5x^ satisfies constraint equation, we have TZ^5x^ = 0. 
We now consider which condition 6x^ must be realized in order to satisfy a constraint equations. 
We can derive the usual relativistic equations of motion only under the condition ( 12 ip . For 
nonholonomic systems a definition of the variations was suggested by Chetaev [HI [20]. The 
variations (5x^ are defined by the condition: 

where 

/ = r]^,uV + c\ (23) 

Using (1211) and fl22|) . we have the functions 7^^ as linear combinations of Of /du^, i.e. 

^^ - ^9^' 

where A is a Lagrange multiplier. We note that substitution of f l23|) into f l22|) gives 

ri^uu^dx" = 0. 
Equations ( 12U]) and ( 12^ give the variational equation 



^ - oD'-^^ir, x,p)-X^]6x^ = 0. (24) 



dr - r - y ^ ^r-j g^^ 

This variational equation is equivalent to the fractional equations of motion 

^=oD2-a^M(r,x,p) + A^ (/i = l,2,3,4). (25) 

We cannot use constraint equation for the function / in variational equation before the partial 
derivative on u'^ is taken. 

Substitution of ( 123|) into ( 125|) gives the equations of motion 

^ = o/^'""-F^(r, X, p) + 2Am^ u^u^ + c^ = 0, (26) 

dr 
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where p'^ = rriQU^ and u'^ = dx'^/dr. The system of equations (l26l) is a closed system of five 
equations in the same number of unknowns x^ and A. Using these equations, we can find the 
multipher A as a function A = A(r, x,p). Substituting this function in (125|) . we get the equations 
for coordinates x^. It allows us to represent the fractional equations of motion for relativistic 
systems with nonholonomic constraint as fractional equations for holonomic systems. 
Differentiating of constraint ([5]) with respect to r, we obtain 

^i^'^-J-P +^o-^c =0. (27) 

Substituting fl27j) in fl26l) with rriou'^ = p^ and dmo/dr = 0, we get 

r7^,p%D^"J-^(r,a;,p) + 2Xr]^,u^p' = 0. 

Using the constraint equation rj^yU^u'^ = —(? and the four-momentum p^ = tuqU^, we obtain 
the Lagrange multiplier 

^ = ^{v,uu\dI-":f^{t,x,p)). 

Therefore the reaction four-force TZ'^ of the nonholonomic constraint is 
As a result, we have the fractional equation 



^=oDl'-:F^{r,x,p) + - 
dr c 



oDl-''J^^{r,x,p) + -Xiu^oDl-'^r). (28) 



These equations define a holonomic system subjected to the sum of four-forces qDI'^^F^ + TZ^ . 
If initial dates satisfy constraint equation (j4]), then the solution of Eq. (!28|) describes a fractional 
dynamics of the relativistic point particle as a holonomic system. 
As a result, we prove the following statement. 

Proposition. Fractional equations for the relativistic particle subjected to a non-potential 
four-force T^ , which have the form 

^ = ±pM, ^=^Dl-'^jr,^T,x,p), r/,.pV + mgc^ = 0. (29) 

or mo dr 



with dmo/dT = 0, are equivalent to the equations 



dr iriQ dr 

where 



p^ -^=oD2-.^/^(r,a;,p) + 7^^(r,x,p), (30) 



1 



n^{T,x,p) = -^p^ {p,,Dl-"Tn, (31) 

and the initial dates satisfy constraint condition Op. 



The solution of Eq. (!30|) describes the fractional dynamics of the relativistic particle. 

5 Fractional non-Hamiltonian and dissipative relativis- 
tic systems 

The system is called locally Hamiltonian if the sum of applied forces satisfies the Helmholtz 
conditions [211 [22]. If (x,p) G M. and A^ is a simply connected region, then a locally Hamil- 
tonian system is globally Hamiltonian. A region is simply connected if it is path-connected 
and every path between two points can be continuously transformed into every other. A region 
where any two points can be joined by a path is called path-connected. 
The Helmholtz conditions for fractional equations (l30l) have the form 



aUDl-Tn^ff^^^^ (32) 



dp" dpu 



dx" dx" dx^^ dx^^ 

Substitution of (EI]) into Eqs. ^ and ([22]) gives 

a(oD^"^^) ^ 1 gK(p.oD^"J--)] _ 
dp" rriQC^ dp" 



0. (33) 



(34) 



dx" tUqC^ \ '^ dx" ) 



ouoi-rn 1 .(.^a^lrfl- ,3,) 



dxi^ mlc^ V dxi^ 
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These equations are the Helinholtz conditions [211 [22] for fractional relativistic dynamics. If 
these conditions are satisfied then the fractional dynamics of relativistic particle is Hamiltonian. 
The fractional relativistic particle subjected to a four-force J^^{t,x,p) is non-Hamiltonian if 
the Helmholtz conditions (IMll and ( l35l) are not satisfied [T5] . 
If 



4 



dioD^-^J^f^) dioD^-'^n^')' 



then we have a generalized dissipative system [15]. If Q{x,p) < for all points {x,p) and 
fl{x,p) < for some points {x,p), then the system is a dissipative system. 

Note that a one-dimensional relativistic particle with dissipation is considered in Refs. [131 
[H]. The Lagrangian and Hamiltonian functions for one-dimensional relativistic particles with 
linear dissipation are suggested. In general, non-Hamiltonian and dissipative n-dimensional 
systems with n > 1 cannot be described by Hamiltonian or Lagrangian since the Helmholtz's 
conditions for these systems are not satisfied [T5] . 

In fractional relativistic dynamics the principle of stationary action for particle subjected 
to non-potential forces J-'^(r, a;,p) can be used if the Helmholtz conditions ( 1M1) and ( l35l) are 
satisfied. The Hamilton's principle and the principle of stationary action are equivalent only for 
special forms of the four- force J-'^{t, x, p). We note that the Hamilton's principle is described by 
nonholonomic variational equation [231 HM [25l [26l [271 [28] . It allows us to use this principle to 
obtain fractional equations of motion for non-Hamiltonian and dissipative systems. The princi- 
ple of stationary action is defined by holonomic variational equation. Therefore the principle of 
stationary action cannot be to derive fractional equations of motion in the general case. In gen- 
eral, the Hamilton's principle and nonholonomic variational equations can be used to describe 
fractional dynamics of relativistic systems. We note that the fractional equations of motion 
which follow from the d'Alembert-Lagrange principle are not equivalent to the fractional equa- 
tions which follow from the principle of stationary action. In Refs. [23l [2ll [251 [221 [30], authors 
give proofs that the solutions to the equations of motion which follow from the d'Alembert- 
Lagrange principle and the Hamilton's principle do not in general satisfy the equations which 
follow from the action principle with nonholonomic constraints. The variational Sedov's equa- 
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tion [261 EH EH] (see also [3ll |32]) can be used in fractional relativistic dynamics instead of 
the principle of stationary action. We note that relativistic models of continuous media with 
dissipation are considered in [281 122] • 

Let us consider the four-vector oD^'^^F^ as the sum 

0^2-" J-^ = G^ + n^, (36) 

where (G^m^) = 0, and [Wu^) ^ 0. Substitution of ([3SD into (EUD of the form ^ gives 

dr & 

The four-force G^ is usually called [17] a real mechanical force, which satisfies the orthogonal 
condition m^C = 0. The four-vector 11'^ describes the energy-momentum exchange between 
the point particle and medium. The components of 11'^ are 

n^ = (7n,(7/c)$), 

where 11 and $ are momentum and energy, which are transmitted by convection per unit time. 
For the heat transfer, three-momentum b'p and energy bQ transmitted per time dr are defined 
by the formulas (5p = Ilfit, and bQ = ^dt. The components of 6Q^ are 

6Q>' = Wdr = (6p, -6Q) = hUdr, -<l>dr), 
c c 

where 6Q'^ is a four-vector of the heat energy-momentum, which is transmitted per time dr. 
Note that the value — m^II'^ = — 7^((n,'u) — $), is a velocity of the convective transmission of 
incoming energy in the rest reference frame. The four-vectors C^ and 11^ allow us to describe 
non-Hamiltonian and dissipative processes in fractional relativistic mechanics. 

6 Conclusion 

We formulate fractional dynamics of relativistic point particle as mechanics of the systems with 
nonholonomic constraint in the four-dimensional pseudo-Euclidean space-time. We consider 
fractional dynamics of relativistic particles subjected to four-forces that can be non-potential. 
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The conditions on the four-forces that allow us to consider fractional dynamics of relativis- 
tic particles subjected to non-potential forces as Hamiltonian dynamics are suggested. We 
prove that the nonholonomic constraint, which represents relativistic invariance, and the non- 
potential four-force can be compensated such that the fractional dynamics is Hamiltonian (and 
non-dissipative) . 

Let us note some possible extensions of the fractional relativistic dynamics. 

1. Nonholonomic constraints with power-law memory [33], which are described by fractional 
equations, can be considered in relativistic mechanics by using fractional derivatives [2] 
with respect to proper time. 

2. The suggested fractional relativistic dynamics can be used to generalize quantum theory 
of non-Hamiltonian and dissipative systems [15] . 



3. In the framework of the fractional relativistic dynamics it is possible to consider a rel- 
ativistic generalization of the fractional variational problems [Ml ES] in Lagrangian and 
Hamiltonain form [36l l37l [38l l39l |40] . Note that nonholonomic variational equations must 
be used since the fractional equations which follow from the d'Alembert-Lagrange princi- 
ple (and the Hamilton's principle) do not in general equivalent the equations which follow 
from the action principle with nonholonomic constraints [23l l24l l25l [29l [30] . 

The study of plasma systems containing ensembles of particles (dust) is a rapidly developing 
field of complex systems research. One of the general features of complex plasma systems is 
the presence of non-potential interaction forces between the dust particles due to the dynamic 
interaction between the dust particles and the plasma (for example, see [H] |l2l US] and refer- 
ences therein). In general, these systems cannot be described as Hamiltonian, since the energy 
is not conserved because of the openness of the systems due to plasma-particle interaction. We 
hope that fractional dynamics of relativistic particle subjected to non-potential forces can be 
used to describe relativistic complex plasma systems. 
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